source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA.R')
plot(f.test)
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA.R')
plot(f.test)
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA.R')
plot(f.test)
x = rlnorm(10^4, mean.vec[i]  , 1)
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA.R')
plot(f.test)
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA.R')
x = rlnorm(10^4, 1,1)
a = data.frame(income = x, group = "A")
x = rlnorm(10^4, mean.vec[i]  , 1)
b= data.table(income = x, group = "B")
d = rbind(a, b)
library(data.table)
x = rlnorm(10^4, 1,1)
a = data.frame(income = x, group = "A")
x = rlnorm(10^4, mean.vec[i]  , 1)
b= data.table(income = x, group = "B")
d = rbind(a, b)
mean = d[, mean(x), by = group]
m = d[, mean(x), by = group]
g.inner = ineq(m$V1)
m = d[, ineq(x), by = group]
ineq(m$V1)
m = d[, ineq(x), by = group]
g.inner = d[, ineq(x), by = group]
g.inner = mean(g.inner$V1)
g.outer = d[, mean(x), by = group]
g.outer = ineq(g.outer$V1)
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA.R')
g.ratio
plot(f.test, g.ratio)
plot(g.ratio)
x = rlnorm(10^4, 1,1)
a = data.frame(income = x, group = "A")
x = rlnorm(10^4, mean.vec[i]  , 1)
b= data.table(income = x, group = "B")
d = rbind(a, b)
anova  = aov(d$income ~ d$group, d)
f = summary(anova)[[1]][[4]]
f.test[i] = f[1]
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA.R')
d = data.table(d)
g.inner = d[, ineq(x), by = group]
i = 2
x = rlnorm(10^4, 1,1)
a = data.table(income = x, group = "A")
x = rlnorm(10^4, mean.vec[i]  , 1)
b= data.table(income = x, group = "B")
d = rbind(a, b)
d = data.table(d)
anova  = aov(d$income ~ d$group, d)
f = summary(anova)[[1]][[4]]
f.test[i] = f[1]
g.inner = d[, ineq(x), by = group]
View(g.inner)
g.inner = mean(g.inner$V1)
g.outer = d[, mean(x), by = group]
x = rlnorm(10^4, mean.vec[i]  , 1)
mean(x)
x = rlnorm(10^4, 1,1)
mean(x)
g.outer = d[, mean(x), by = group]
x = rlnorm(10^4, 1,1)
a = data.table(income = x, group = "A")
x = rlnorm(10^4, mean.vec[i]  , 1)
b= data.table(income = x, group = "B")
d = rbind(a, b)
d = data.table(d)
g.inner = d[, ineq(x), by = group]
g.inner = mean(g.inner$V1)
g.outer = d[, mean(x), by = group]
g.outer = ineq(g.outer$V1)
View(d)
g.inner = d[, ineq(income), by = group]
View(mean)
g.inner = d[, ineq(income), by = group]
g.inner = mean(g.inner$V1)
g.outer = d[, mean(income), by = group]
g.outer = ineq(g.outer$V1)
g.ratio[i] = g.outer/g.inner
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA.R')
plot(f.test, g.ratio, log = "x")
plot(f.test, g.ratio, log = "xy")
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA.R')
plot(f.test, g.ratio)
lm(log(f.test) ~ log(g.ratio))
r = lm(log(f.test) ~ log(g.ratio))
summary(r)
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA.R')
y <- rnorm(n = 100)
x <- matrix(data = rnorm(1000), nrow = 100)
a <- 1:10
fast.anova(x = x, y = y, assign = a)
install.packages("hit")
library(hit)
y <- rnorm(n = 100)
x <- matrix(data = rnorm(1000), nrow = 100)
a <- 1:10
fast.anova(x = x, y = y, assign = a)
fast.anova(x = x, y = y)
View(x)
y <- rnorm(n = 100)
x <- matrix(data = rnorm(1000), nrow = 100)
a <- 1:10
fast.anova(x = x, y = y, assign = a)
library(matrixStats)
library(hit)
colVars(x)
x = rnorm(10^5, 1,1)
sum( (x - mean(x) )^2 )
var(x) * (length(x) - 1)
library(microbenchmark)
microbenchmark(
sum( (x - mean(x) )^2 ),
var(x) * (length(x) - 1)
)
source('~/.active-rstudio-document')
microbenchmark(
sum( (x - mean(x) )^2 ),
var(x) * (length(x) - 1)
)
microbenchmark(
sum( (x - sum(x)/length(x) )^2 ),
var(x) * (length(x) - 1)
)
microbenchmark(
sum( (x - sum(x)/length(x) )^2 ),
var(x) * (length(x) - 1)
)
microbenchmark(
sum( (x - sum(x)/length(x) )^2 ),
var(x) * (length(x) - 1)
)
microbenchmark(
sum( (x - sum(x)/length(x) )^2 ),
var(x) * (length(x) - 1)
)
microbenchmark(
sum( (x - sum(x)/length(x) )^2 ),
var(x) * (length(x) - 1),
times = 500
)
microbenchmark(
sum( (x - sum(x)/length(x) )^2 ),
var(x) * (length(x) - 1),
times = 500
)
sum( (x - sum(x)/length(x) )^2 )
var(x) * (length(x) - 1)
n = 10^5
x = rnorm(n, 1,1)
y = seq(1,n, length.out = 10 )
y
y = rep(1:n.group, each = n/n.group)
n = 10^5
n.group = 10
x = rnorm(n, 1,1)
y = rep(1:n.group, each = n/n.group)
library(data.table)
d = data.table(x, y)
d = data.table(x, y = as.factor(y))
t1 = aov(x ~ y, d)
summary(t1)
ss = function(x){  sum( (x - sum(x)/length(x) )^2 )  }
ss(x)
test = d[, ss(x), b = y]
View(test)
test = d[, ss(x), b = y]
View(test)
summary(t1)
n = 10^5
n.group = 10
x = rnorm(n, 1,1)
y = rep(1:n.group, each = n/n.group)
d = data.table(x, y = as.factor(y))
t1 = aov(x ~ y, d)
summary(t1)
sum(d$x)
n = 10^5
n.group = 10
x = rnorm(n, 1,1)
y = rep(1:n.group, each = n/n.group)
d = data.table(x, y = as.factor(y))
t1 = aov(x ~ y, d)
summary(t1)
ss = function(x){  sum( (x - sum(x)/length(x) )^2 )  }
test = d[, ss(x), b = y]
sum(d$x)
n = 10^5
n.group = 10
x = rnorm(n, 1, 2)
y = rep(1:n.group, each = n/n.group)
d = data.table(x, y = as.factor(y))
t1 = aov(x ~ y, d)
summary(t1)
ss = function(x){  sum( (x - sum(x)/length(x) )^2 )  }
test = d[, ss(x), b = y]
sum(d$x)
n = 10^5
n.group = 10
x = rnorm(n, 1, 2)
y = rep(1:n.group, each = n/n.group)
d = data.table(x, y = as.factor(y))
t1 = aov(x ~ y, d)
summary(t1)
ss = function(x){  sum( (x - sum(x)/length(x) )^2 )  }
test = d[, ss(x), b = y]
sum(d$x)
View(test)
sum(d$y)
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA Function.R')
View(d)
sum(d$x)
sum(test$y)
sum(test$V1)
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA Function.R')
summary(t1)
sum(test$V1)/10^4
sum(test$V1)/99990
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA Function.R')
summary(t1)
sum(test$V1)/99990
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA Function.R')
test = d[, ss(x), b = y]
summary(t1)
sum(test$V1)/99990
test = d[, mean(x), b = y]
test2 = ss(test$V1)
View(d)
View(test)
group.size = n/n.group
test2 = ss(test$V1)*group.size
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA Function.R')
test2 = ss(test$V1)*group.size
summary(t1)
test = d[, mean(x), b = y]$V1
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA Function.R')
ss.between  = group.size*ss(d[, mean(x), b = y]$V1)
summary(t1)
ss.total = ss(x)
eta.square = ss.between/ss.total
eta.square = group.size*ss(d[, mean(x), b = y]$V1)/ss(x)
summary(t1)
library(lsr)
etaSquared(t1)
etaSquared(t1) - eta.square
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA Function.R')
etaSquared(t1) - eta.square
etaSquared(t1) - eta.square
etaSquared(t1) - eta.square
etaSquared(t1) - eta.square
eta.r = etaSquared(t1)
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA Function.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA Function.R')
eta.r = etaSquared(aov(x ~ y, d))[1]
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA Function.R')
microbenchmark(
eta.square = group.size*ss(d[, mean(x), b = y]$V1)/ss(x),
eta.r = etaSquared(aov(x ~ y, d))[1],
times = 10)
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA Function.R')
microbenchmark(
eta.square = group.size*ss(d[, mean(x), b = y]$V1)/ss(x),
eta.r = etaSquared(aov(x ~ y, d))[1],
times = 10)
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA Function.R')
microbenchmark(
eta.square = group.size*ss(d[, mean(x), b = y]$V1)/ss(x),
eta.r = etaSquared(aov(x ~ y, d))[1],
times = 4)
n = 10^5
n.group = 500
group.size = n/n.group
x = rnorm(n, 1, 2)
y = rep(1:n.group, each = n/n.group)
d = data.table(x, y = as.factor(y))
#t1 = aov(x ~ y, d)
ss = function(x){  sum( (x - sum(x)/length(x) )^2 )  }
ss.between  = group.size*ss(d[, mean(x), b = y]$V1)
ss.total = ss(x)
eta.square = ss.between/ss.total
microbenchmark(
eta.square = group.size*ss(d[, mean(x), b = y]$V1)/ss(x),
eta.r = etaSquared(aov(x ~ y, d))[1],
times = 1)
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA Function.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA Function.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA Function.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA Function.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA Function.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA Function.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA Function.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA Function.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA Function.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/ANOVA Function.R')
setwd("~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test")
p = read_csv("Parameters.csv")
library(ineq)
library(data.table)
library(NORMT3)
library(blair)
library(microbenchmark)
library(Rcpp)
library(RcppZiggurat)
library(snow)
library(doSNOW)
library(Runuran)
library(magrittr)
library(readr)
p = read_csv("Parameters.csv")
View(p)
i = 1
n.groups = p$n.groups[i]
group.size = round(10^5/n.groups)
n = group.size*n.groups
sd.mu = p$sd.mu[i]
mu =  rlnorm(n.groups, 1, sd.mu)
mu = rep(mu, each = group.size)
m.sigma =  p$m.sigma
m.sigma =  p$m.sigma[i]
sd.sigma = p$sd.sigma[i]
mu =  rlnorm(n.groups, 1, sd.mu)
mu =  urlnorm(n.groups, 1, sd.mu)
microbenchmark(
mu =  urlnorm(n.groups, 1, sd.mu),
mu2 =  rlnorm(n.groups, 1, sd.mu)
)
microbenchmark(
mu =  urlnorm(10^5, 1, sd.mu),
mu2 =  rlnorm(10^5, 1, sd.mu)
)
microbenchmark(
mu =  urlnorm(4, 1, sd.mu),
mu2 =  rlnorm(4, 1, sd.mu)
)
microbenchmark(
mu =  rlnorm(4, 1, sd.mu),
mu2 =  exp(zrnorm(4)*sd.mu +1   )
)
microbenchmark(
mu =  rlnorm(4, 1, sd.mu),
mu2 =  exp(zrnorm(4)*sd.mu +1   )
)
n.groups = p$n.groups[i]
group.size = round(10^5/n.groups)
n = group.size*n.groups
sd.mu = p$sd.mu[i]
mu =  rlnorm(n.groups, 1, sd.mu)
mu = rep(mu, each = group.size)
m.sigma =  p$m.sigma[i]
sd.sigma = p$sd.sigma[i]
sigma = rlnorm(n.groups, m.sigma, sd.sigma)
sigma = rep(sigma, each = group.size)
x = exp(zrnorm(n)*sigma + mu)
group.id = rep(1:n.groups, each = group.size)
d = data.table(x, group.id)
g.inner = d[, fastgini(x), by = group.id]
g.inner = mean(g.inner$V1)
g.outer = d[, mean(x), by = group.id]
g.outer = fastgini(g.outer$V1)*(n.groups)/(n.groups-1)
g.ratio = g.outer/g.inner
g = fastgini(x)
data.table(g.inner, g.outer, g.ratio,  g, n.groups )
eta.square = group.size*ss(d[, mean(x), b = y]$V1)/ss(x)
ss = function(x){  sum( (x - sum(x)/length(x) )^2 )  }
eta.square = group.size*ss(d[, mean(x), b = y]$V1)/ss(x)
eta.square = group.size*ss(d[, mean(x), b = group.id]$V1)/ss(x)
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/Model.R')
plot(good$g.ratio, good$eta.square)
View(p)
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/Generate Parameters.R')
View(good)
plot(good$g.ratio, good$eta.square)
plot(good$g.ratio, good$eta.square, log = "x")
plot(good$g.ratio, good$eta.square, log = "xy")
runTime
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/Generate Parameters.R')
plot(good$g.ratio, good$eta.square)
plot(good$g, good$eta.square, log = "xy")
plot(good$g.outer, good$eta.square, log = "xy")
good = subset(model, abs(model$g -0.5) <= 0.02 ) %>% na.omit
#plot(good)
#write_csv(good, "Parameters.csv")
plot(good$g.outer, good$eta.square, log = "xy")
good = subset(model, abs(model$g -0.5) <= 0.05 ) %>% na.omit
#plot(good)
#write_csv(good, "Parameters.csv")
plot(good$g.outer, good$eta.square, log = "xy")
good = subset(model, abs(model$g -0.5) <= 0.01 ) %>% na.omit
#plot(good)
#write_csv(good, "Parameters.csv")
plot(good$g.outer, good$eta.square, log = "xy")
plot(good$g.ratio, good$eta.square, log = "xy")
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/Generate Parameters.R')
good = subset(model, abs(model$g -0.5) <= 0.025 ) %>% na.omit
plot(good$g.ratio, good$eta.square, log = "xy")
good = subset(model, abs(model$g -0.5) <= 0.05 ) %>% na.omit
plot(good$g.ratio, good$eta.square, log = "xy")
plot(good)
write_csv(good, "Parameters.csv")
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Metric Test/Cohen D.R')
plot(g.ratio, d)
source('~/Desktop/Empirical Research/Plots/Group Income Effect/Eta_plot.R')
exp(coef(r))
source('~/Desktop/Empirical Research/Plots/Group Income Effect/Eta_plot.R')
source('~/Desktop/Empirical Research/Plots/Group Income Effect/Eta_plot.R')
source('~/Desktop/Empirical Research/Plots/Group Income Effect/Eta_plot.R')
source('~/Desktop/Empirical Research/Plots/Group Income Effect/Model Benchmark.R')
source('~/Desktop/Empirical Research/Plots/Group Income Effect/Model Benchmark.R')
source('~/Desktop/Empirical Research/Plots/Group Income Effect/Model Benchmark.R')
source('~/Desktop/Empirical Research/Plots/Group Income Effect/Model Benchmark.R')
View(m)
source('~/Desktop/Empirical Research/Plots/Group Income Effect/Model Benchmark.R')
View(m)
setwd("/home/blair/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Model Bootstrap 2")
m = read_csv("Firm Bootstrap Results.csv")
m = data.table(m)
View(m)
View(h)
h = read_csv("Level Bootstrap Results.csv")
h = data.table(h)
h.mean = h[ , as.list(get.summary(mean)), by = level]
View(h.mean)
h.gini = h[ , as.list(get.summary(gini)), by = level]
View(h.gini)
View(h)
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Model Bootstrap 2/Bootstrap Script Big Foreach RCPP.R')
source('~/Desktop/Empirical Research/Plots/Group Income Effect/Model Benchmark.R')
source('~/Desktop/Empirical Research/Plots/Group Income Effect/Boxplots II.R')
source('~/Desktop/Empirical Research/Plots/Group Income Effect/Boxplots II.R')
View(final)
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Occupation/BLS/BLS Analysis III.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Occupation/BLS/FUNCTIONMASTER SCRIPT.R')
dir = dirname(rstudioapi::getActiveDocumentContext()$path)
dir = dirname(rstudioapi::getActiveDocumentContext()$path)
workspace = gsub("Functions", "BLS Data", dir)
setwd(workspace)
O = read_csv("Occupations.csv")
dir  = dirname(rstudioapi::getActiveDocumentContext()$path)
workspace = paste(dir, "/Functions", sep = "")
setwd(workspace)
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Occupation/BLS/Functions/burr.R')
stopCluster(cl)
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Occupation/BLS/FUNCTIONMASTER SCRIPT.R')
dir  = dirname(rstudioapi::getActiveDocumentContext()$path)
workspace = paste(dir, "/Functions", sep = "")
workspace
dir
dir  = dirname(rstudioapi::getActiveDocumentContext()$path)
workspace = paste(dir, "/Functions", sep = "")
setwd(workspace)
list.files(".R")
list.files()
list.files(pattern = ".R")
dir  = dirname(rstudioapi::getActiveDocumentContext()$path)
workspace = paste(dir, "/Functions", sep = "")
setwd(workspace)
list.files(pattern = ".R")
list.files(pattern = ".R")
list.files(pattern = ".R ")
list.files(pattern = ".R")
list.files(pattern = ".csv")
list.files(pattern = ".R~")
list.files(pattern = "\\.R$")
f = list.files(pattern = "\\.R$")
source(f[1])
dir  = dirname(rstudioapi::getActiveDocumentContext()$path)
dir
dir = dirname(rstudioapi::getActiveDocumentContext()$path)
workspace = paste(dir, "/Functions", sep = "")
setwd(workspace)
files = list.files()
dir  = dirname(rstudioapi::getActiveDocumentContext()$path)
files = list.files()
files
files = list.files(pattern = "\\.R")
list.files(pattern = "\\.R")
list.files(pattern = "\\.R$")
files = list.files(pattern = "\\.R$")
files = files[files != "FUNCTIONMASTER SCRIPT.R"]
files
source(file[3])
source(files[3])
stopCluster()
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Occupation/BLS/Functions/FUNCTIONMASTER SCRIPT.R')
library(R.utils)
dir  = dirname(rstudioapi::getActiveDocumentContext()$path)
files = list.files(pattern = "\\.R$")
files = files[files != "FUNCTION MASTER SCRIPT.R"]
files = list.files(pattern = "\\.R$")
files = files[files != "FUNCTION MASTER SCRIPT.R"]
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Aggregate Studies/Mueller/Model/Mueller Analysis.R')
files = files[files != "Mueller.csv"]
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Aggregate Studies/Mueller/Model/Mueller Analysis.R')
